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Abstract 

The long time behavior of an absorbed Markov process is well described by the 
limiting distribution of the process conditioned to not be killed when it is observed. 
Our aim is to give an approximation's method of this limit, when the process is a 1- 
dimensional Ito diffusion whose drift is allowed to explode at the boundary. In a first 
step, we show how to restrict the study to the case of a diffusion with values in a bounded 
interval and whose drift is bounded. In a second step, we show an approximation 
method of the limiting conditional distribution of such diffusions, based on a Fleming- 
Viot type interacting particle system. We end the paper with two numerical applications 
: to the logistic Feller diffusion and to the Wright-Fisher diffusion with values in ]0,1[ 
conditioned to be killed at 0. 
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1 Introduction 

Let (X t ) be a killed Markov process with law P, taking its values in E U {d}, where d is a 
cemetery point. We denote by tq = inf {t > 0, X t = d} the killing time of (X t ). A probability 
measure v on E is called a quasi-stationary distribution (QSD) if, for all t > 0, the 
distribution of the process X, initially distributed with respect to v and conditioned to be 
not killed before time t, is still v at time t, that is F u (X t e A\tq > t) = u(A) for every 
A C E and t > 0. Without loss of generality, we suppose that d is an absorbing point, so 
that {r d > t} = {X t ^ d}. 

Let /i be a probability measure on E. If it exists and provided it is a probability, the 
limiting conditional distribution 

lim F x (X t G \X t ±d) 

t— >+oo 
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is called the Yaglom limit for /x, from the Russian Mathematician A.M. Yaglom. He showed 



in 



27| that the limiting conditional distribution of the number of descendants in the 



n 



generation of a Galton- Watson process always exists in the subcritical case. 

The existence or uniqueness of such invariant conditional distributions have been proved 
in a host of contexts. When E is finite, it is proved in [t| that there exists a unique QSD 
v and that the Yaglom limit converges to v independently of the initial distribution. In jij, 
the case of a birth and death process on N is studied. For this process, the set of QSDs is 
either empty, or a singleton, or a continuum indexed by a real parameter and given by an 
explicit recursive formula. This is an exception : most of the known results on QSDs are 



related with existence or uniqueness problems. In [111 ], the existence of a quasi-stationary 
distribution for a continuous time Markov chain on N killed at is proved under conditions 
on moments of the killing time, using an original renewal dynamical approach. In the 
case of 1-dimensional diffusion on [0, + oo[ with C 1 drift and killed at are studied, with the 
assumption that +oo is a natural boundary. The dependence between the initial measure 



and the Yaglom limit is explored in [19( (for a Brownian motion with constant drift killed at 



0) and 18j (for the Orstein-Uhlenbeck process killed at 0). In 26|], the case of 1-dimensional 
diffusions with general killing on the interior of a given interval is investigated. In Q|, the 
authors study the existence and uniqueness of the QSD for 1-dimensional diffusions killed at 
and whose drift is allowed to explode at the boundary, which is the case under study in the 



present paper. See [23| for a regularly updated extensive bibliography on QSD. 

In this paper we are concerned with 1-dimensional Ito diffusions with values in ]0, + 
oo[U{<9} killed at and defined by the stochastic differential equation 

dX t = dB t - q{X t )dt, X = x > 0, 

where B is a standard 1-dimensional Brownian motion and q G C 1 (]0, + oo[). In [3[, the 
Yaglom limit of this process is studied and the authors give some conditions on the drift q, 
which are sufficient for the existence and the uniqueness of the QSD. In particular, they allow 
the drift to explode at the origin. As explained in the paper, this diffusion is closely related 
with some Markov mortality models. Such applications need the computation of the process 
QSD, but the tools used in [3( are based on spectral theory's arguments and don't allow us 
to get explicit values. Our aim is to give an easily simulable approximation's method of this 
QSD. 

The problem of QSD's approximation has been already explored in 0], 12] when E is 



a bounded open set of M d and X is a Brownian motion killed at the boundary of E. The 
authors proved an approximation's method exposed in [l|, which is based on a Fleming- 
Viot type system of interacting particles whose number is going to infinity. In 10J, it is 
proved that this method works well for a continuous time Markov chain in a countable 
state space under suitable assumption on the transition's rates (moreover, the existence of 
a QSD is a consequence of the approximation's method). New difficulties arise from our 
case with unbounded drift. For instance, the interacting particle process introduced in [l[ 
isn't necessarily well defined. To avoid this difficulty, we begin by proving that one can 
approximate our QSD by the QSD's of diffusions with bounded drifts. 
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Let us denote by F e the law of a diffusion with values in ]e,l/e[, defined by the stochastic 



differential equation dX t = dB t — q(X t )dt and killed when it hits e or f/e. In 22], it is proved 



that the Yaglom limit associated with P e exists and is its unique QSD. We will denote it by 
z/ e . In the first part of this paper, we give some conditions on q G C 1 (]0, + oo[) for the family 
(^e)o<e<i/2 to be tight and to converge, when e — > 0, to a QSD for the law P°. We point 
out the fact that this result remains valid in the case of an unbounded drift diffusion with 
values in a bounded interval. In a second part, we prove an approximation method for each 
probability measure u e , based on the interacting process introduced in Fix e > and let 
us describe the interacting particle process of size N > 2: each particle moves independently 
in ]e,l/e[, each one with law P e until one of them hits the boundary. At this time, the killed 
particle jumps on the position of an other particle, chosen uniformly between the N — 1 
remaining one. Then the particles evolve independently, until one of them is killed and so on 
(see Figured]). One has to prove that the particles don't degenerate at the boundary. In 0], 




Figure 1: The interacting particle system (X l ,X 2S] 



the authors prove a non- degeneracy result with arguments based on a construction of the d- 
dimensional Brownian motion due to Ito, where d > 2. It seems that this tool can't be easily 
generalized to other diffusions. To prove such results under our settings, we build an original 
coupling between the interacting particle process and an independent particle system of the 
same size. This coupling is valid for all drifted Brownian motions with continuous bounded 
drift, killed at the boundary of a bounded interval of M. + . It will be used in each step of the 
proof. 

We conclude the paper by two numerical applications. At first, we treat the case of the 
logistic Feller diffusion introduced in [17| and studied in [3] with values in ]0, + oof, driven 
by the stochastic differential equation 

dZ t = y^tdB t + (rZ t - cZf)dt, Z = z > 0, 
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where B is a 1- dimensional Brownian motion and r, c are two positive constants. Clearly, 
is an absorbing state for this diffusion. In a second time we study in detail the case of the 
Wright-Fisher diffusion on ]0,1[ conditioned to be killed at I 
values in ]0,1[, and is defined by 



see 



13]). This diffusion takes 



dZ t = yJZ t {\ - ZtjdB t + (1 - Z t )dt, Z = z e]0,l[, 
where B is a 1-dimensional Brownian motion. This diffusion is absorbed at 1. 



2 From unbounded drift to bounded drift 

Let P° be the law of a diffusion process taking its values in ]0, + oo[U{<9}, killed when it hits 
and defined by the stochastic differential equation (SDE) 

dX t = dB t - q(X t )dt, X = x > 0, 

where B is a 1-dimensional Brownian motion. The drift q is taken in the set of real valued 
continuously differentiable functions C 1 (]0,+oo[). We denote by L° the infinitesimal generator 
associated with P°. 

We define, Vx e]0, + oo[, 



Q{x) = J q(y)dy, 
dfx(x) = e~ 2Q{x) dx, 



and 



W{x) = q(x) 2 — q'(x). 

For all e e]0,l/2[, we define F e x as the law of the diffusion taking its values in ]e,l/e[, 
defined by the SDE 

dX t = dB t - q(X t )dt, X = x e]e,l/e[ 

and killed when it hits the boundary {e,l/e}. Let L e be the infinitesimal generator 1/2 A — qV 
with the Dirichlet boundary condition on {e,l/e}. — L e has a simple real eigenvalue X e (see 



22l . Theorem KR]) at the bottom of its spectrum. The corresponding eigenfunction r] e is 

e 

rj t (x) 2 dfj,(x) 



positive and belongs to C 2 ([e,l/e]). We choose it so that 

•1A 



Let us recall some results of |22|: 

Theorem (Pinsky (1985)) The Yaglom limit associated with P e exists for all initial dis- 
tributions 5 X , x e]e,l/e[, and doesn't depend on x. This limit is a QSD, which we denote by 
u e . Furthermore, we have 

rj e {x)diJi{x) 



dv e (x) 



fe /€ Ve(x)dn(x) 
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In fact, u e is the unique QSD of the process, as proved in Lemma [T9l below, but we won't use 
it in this section. The aim of this section is to study the asymptotic behaviour of (u e ) when 
e goes to 0. 

From now, .Mi(]0, + oo[) denotes the space of probability measures on ]0, + oof equipped 
with the weak topology. The following hypotheses arise naturally in the proof of Theorem [TJ, 
which is based on a compactness-uniqueness method. 

Hypothesis 1 (HI) W is bounded below by —C, where C is a positive constant. Moreover, 
W(x) — > +00 when x — * 00. 



Theorem 1 Assume that hypotheses (HI) and (H2 or H3) are satisfied. Then 

v e — >u E Mi(]0, + oo[), 

where v is a QSD for P°, which is equal to the Yaglom limit lim^ +00 P°(Xi G .\t < tq), 
Vx G]0, + oo[. 

Remark 1 The hypotheses (HI) and (H2 or H3) are the assumptions that are made in 0] 
to prove the existence of the Yaglom limit. 

Remark 2 If a process satisfies the hypotheses of Theorem [TJ then it is killed in finite time 
a.s or it is never killed a.s. Indeed, assume that the process can be killed in finite time with a 
positive probability. Then f? (f* e~^^dy) dx < +00 (see [HI) and J^°° e Q ^dx = +00 
(as a consequence of (HI) and (H2 or H3)). But this two conditions are fulfilled if and only 
if the process is killed in finite time almost surely (see 0, Theorem 3.2 p.450]). 

Remark 3 The existence of a QSD for P° can be seen as a consequence of Theorem [TJ The 
existence of the Yaglom limit is proved in [3j, Theorem 5.2]. 

Remark 4 In Part 12.31 we give the counterpart of Theorem [TJ for diffusions with values in a 
bounded interval. 

The end of the section is devoted to the proof of Theorem [TJ 



Hypothesis 2 (H2) 




Hypothesis 3 (H3) 




e 2Q ^ x) dx < +00 and 
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2.1 Tightness of the family (^ e )o< e <i/2 
This part is devoted to the proof of the following result, 

Proposition 2 Assume that the hypotheses (HI) and (H2 or H3) are satisfied. Then the 
family (v e )o<e<i/2 is tight. Moreover, every limit point is absolutely continuous with respect 
to the Lebesgue measure. 

We know that L e r] e = — \ e r] t , rj e e C 2 ([e,l/e]), and r] t satisfies the differential equation 



l/2<(x) - q(xX{x) 




0. 




1. 



We have 



v e (x)W(x) - v'l{x) = 2X e v e (x), 



with the boundary conditions 



v e (e) = v e {l/e) = 0. 



(3) 



Lemma 3 Assume that the hypothesis (HI) is fulfilled. Then (f e )o<e<i/2 is uniformly bounded 
above and the family (v^(x)dx) 0<t< i/2 is tight. 



Proof of Lemma\^ : From the differential equation satisfied by v e , we have 



v t (x) 2 W(x) -vl(x)v t (x) = 2X e v t (x) 2 . 



Integrating by parts and looking at the boundary conditions ([3]), 




where v e is normalized in L 2 



(dx). That implies 




The eigenvalue X e of — L e is given by (see for instance [28|, chapter XI, part 8]) 




(4) 



6 



where (]e,l/e[) is the vector space of infinitely differentiable functions with compact sup- 
port in ]e,l/e[ and (/,<?) M = / +o ° f(u)g(u)dn(u). We deduce from it that X e increases with e 
and is uniformly bounded above by X 1 / 2 . 
We have then 

fl/e rl/e 

0< / v' e (x) 2 dx + / v e (x) 2 (W(x) + C + l)dx<2\ 1/2 + C+l. (5) 

Looking at the boundary conditions (j3J), we obtain, for all x e]e,l/e[, 

/■iA 

w e 2 (x) = -2 / v' € {y)v e {y)dy 

J X 

2 



< : / v' c (y)v c {y)y/W(y) + C + ldy. 

nun [a . 1/e[ VW 7 + C + 1 A 

Then, applying the Cauchy-Schwarz inequality to the right term above, 



rfU) < " W l 1 \>(y)*dyj [ 1/ \ e (^(\V(! 1 ) + C+ I W/y. 

min M/e[ Vlf + C + l 

From (j3J), the integral product is bounded by 2Ai/2 + C + 1, thus > 0, independent from 
e, such that 

v 2 e (x) < A 



mm [xA/e[ y/W + C + 1 

< ; (6) 

min[ X)+00 [ + C + 1 

where W(x) + C + 1 > 1 for all x e]0, + oo[, thanks to Hypothesis (HI). That implies the 
first part of Lemma [31 

Let us prove that the family (v 2 dx) 0<€< i/ 2 is tight. Fix 5 > 0. We have to find a compact 
subset Kg in ]0, + oo[ such that 



/ 

J]o,- 



v 2 (x)dx < 5, (7) 

-oo[\K s 



for all e e]0,l/2[. Thanks to ®, we have v 2 < A, then 

.S/(2A) 

v 2 < 5/2. 



From the second part of Hypothesis (HI), 3M 5 > such that W(x) + C+1 > 2(2\ + C+l)/5 
for all x > Ms. That implies 

v e {xfdx < < x ) 2 5 %2\Vc + ^ dX 

< 6/2, 

where the last inequality is due to (|3j). Finally, the compact set K$ = [6/ (2 A), Ms] satisfies 
©. □ 



M6 



7 



Lemma 4 Assume that (HI) is satisfied. Then f 1 ^ 6 r] e {y)d^{y) is uniformly bounded below 
by a constant B > 0. 

Proof of Lemma \4\ ■ Assume that f r) e (y)d/j l (y) isn't uniformly bounded below : one can 
find a sub-sequence f^ €k v €k (y)e~ Q( - y ^dy = f^ tk r] ek (y)dfi(y), where — > 0, which tends to 

0. From Lemma [3], (u e )o<e<i/2 is uniformly bounded, so that / e v ek {y) 2 e~ Q ^dy — > 0. The 
family (v e (x) 2 dx) being tight, one can find (after extracting a sub-sequence) a positive map 
m such that, for all continuous and bounded : R + — > R, 

l/e* /-+oo 

Ve k {y) 2( t>{y) d y / ™{y)<l>{y)dy. (8) 

Indeed, (i^ ) being uniformly bounded, all limit measure is absolutely continuous with respect 
to the Lebesgue measure. In particular, 

/■l/ejt /-+00 

/ f efc (y) 2 min(e~ Q(?/) ,l)<i?/ — ► / 771(2/) min (e~ Q( - y \l)dy, 

Je k JO 



then 



P+OO 

/ m(y) min (e~ Q( * y \l)dy = 0. 



But min (e ^,1) is continuous and positive on K + , so that m vanishes almost every where. 
Finally, by the convergence property jS]) applied to equal to 1 almost everywhere, we have 



v f dx — ► 0, 



1/6* 

2 

(.' 



what is absurd. Thus, one can define B = inf £ J e /e rj e (y)dfi(y)/A > 0. □ 

Lemma 5 Assume that (HI) and (H2) are satisfied. Then the family (r) e (x)d/j,(x))o <e< i/2 is 
tight. 

Proof of Lemma [H : By (jSJ) , we have 

l/ \ 2 £ (y)(W(y) + C + l)dti(y) = v 2 e (y)(W(y) + C + l)dy 

< Xl/2 + C + l. 

For all 5,M > 0, using Cauchy-Schwarz inequality, we get on one hand 

<5 / rS \ 5 / r& 

2 



ri,(y)Mv) < (/ n,(y)\w{ v ) + c + i)Mv) I 1 / w(v) + c + i My) > l " 1 



5 (v + c + i)M/ o H , (y) + C+1 ^)) . do) 



On the other hand, 



»7 e (y)^(2/) < / Ve(y)dfi(y) / (n) 

M \J M J \JM J 

/ P+OQ \ § 

< (J d^y)\ . (12) 

Thanks to (H2), both terms are going to uniformly in e, when 5 and M tend respectively 
to and +oo. As a consequence, the family (r] e (x)dfi(x)) 0<e< i/2 is tight. □ 

Lemma 6 Assume that (HI) and (H3) hold. Then the family (T) e { x )d[i(x))o<e<i/2 is tight. 

Proof of Lemma® : The first part of the hypothesis (H3) is the same as (H2)'s one, then 

+oo 

Ve(y)dfi(y) -> 

At 

when M goes to infinity, uniformly in e. 

Moreover, there exists a constant K > such that, for any x G]0,1] and any e g]0,1/2], 

r] € (x) < Kxe Q(x) . 

This is a consequence of 0, Proposition 4.3] whose proof is still available under our settings. 
This inequality allows us to conclude the proof of Lemma [6j □ 

Thanks to equality (j2J) and Lemmas HI [5] and [61 the first part of Proposition [2] is proved. 
Moreover, u e has a density with respect to the Lebesgue measure which is bounded on every 
compact set, uniformly in e > 0. Thus every limit point is absolutely continuous with respect 
to the Lebesgue measure. 

2.2 The limit points of the family (^ e )o< e <i/2 

Proposition 7 Assume that Hypotheses (HI) and (H2 or H3) are fulfilled and let v be a 
probability measure which is the limit of a sub-sequence (v^^n, where — > when k — » oo. 
Then v is a QSD with respect to P°. 

Proof of Proposition [?] : From Proposition [2, the family (^ e )o<e<i/2 is tight. Let v be a 
limit point of the family (^ e )o<e<i/2- There exists a sub-sequence (z^Jfc which converges to u, 
where {ek)km is a decreasing sequence which tends to 0. We already know that v is absolutely 
continuous with respect to the Lebesgue measure. That implies that, for all open intervals 
D =]c,d[c R+, 

v ek (D)^v(D), (13) 
and, for all bounded maps <fi continuous on R + , 

4>(x)dis ek (x) — > / (f)(x)du(x). (14) 
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Let v t (resp. z/ e> t) be the distribution at time t of a diffusion with law P° (resp. P£ ), 
conditioned to be not killed until time £, that is 

^(c/a;) = F° u (uJt e dx\r d > t) 

and 

i/e.tCtte) = KX u t £dx\r d >t) 

The probability measure v e being a QSD for P e , we have u Et t = v e for all t > 0. We want to 
show that v = v t for allt > 0, we have then to prove the following convergence result: 

Vt > 0, =]c,d[C R+, v ekit (D) — h. !/ t (D). (15) 

fc— >+oo 

Indeed, suppose that (|T5|) holds, then on the one hand, v €k {D) = v ekt {D) — > u t (D). In the 
other hand v €k (D) — > v{D). We have then u t (]c,d[) = v(]c,d[), V\c,d[C M + and this conclude 
the proof of Proposition 

Let us prove ffT5l) . By definition, 



v e k ,t{D) 



lH,i/e k[ K k (^eD)du ek (x 

I]e k ,l,e k[ ^9>t)du ek {x) 



The numerator is equal to 

/ G D)du ek (x) = [ ¥° x (u; t G D)dv ek {x) 



+ / [¥^(u H eD)-V p x (u H eD)]du tk (x) 

J]e k d/e k \ 



'Vkdhk 

For all t > 0, the map x t— > P°(cu 4 G D) is continuous and bounded, then, by the convergence 
property ffT4l) . 

/ P°(u; t G D)dv ek (x) -> / P°(a; t G D)dv{x), 

Assume that (H2) is fulfilled, then, similarly to (Q and (jlip . we have, for all bounded 
continuous functions / :]0, + oo[— > R and all M > 0, 

i 

+oo / r+oo \ 2 

f(x)ij ek (x)dfi(x) < [ \f(x)\ 2 dn(x 

M \J M 

and, for all m > 0, 



f(x) Vek (xW(x)<{\ 1/2 + C+lf [ I w ff+c + ^ ti*) 



10 



Replacing f(x) by P°(c^t E D) — P|*(o>t G -D), which is decreasing to when — > oo, and by 
monotone convergence theorem, we have 

/ [F?(t* E D) — F° x (uj t E D)] du ek (x) —> 



and 



/ [P?(<* ED)- F° x (uj t E D)] du th (x)—> 0. 



Finally, the density of v tk being bounded above in every compact set [m,M], uniformly in e*., 
the same argument of monotone convergence gives us 

/ [FZfa ED)- ¥ x (uj t E D)] dv ek (x)—> 0. 

With similar arguments, the same holds under (H3). Finally, we obtain 

/ W*{u t E D)du ek {x) -> f P> t E D)dv{x). 

J]e k ,l/e k [ J]0,+co[ 

Thanks to 0, Lemma 5.3 and Theorem 2.3], the map x i— > V x (tq > t) = P°(^ 4 g]0, + oo[) 
is continuous, and F x (rg > t) — P^fVg > t) is increasing to when k — > oo. Thus the 
denominator can be treated in the same way. □ 

We can now conclude the proof of Theorem [TJ 

Proposition 8 Assume that (HI) and (H2 or H3) hold. The limit measure v in the state- 
ment of Proposition^ is unique. Moreover v is the Yaglom limit associated with P° ; Vx G 
]0, + oo[. 

Proof of Proposition [3 : The proof of Proposition [7] implies that 

K t (r a >t)— >K(ra>t), V* > 0. 
The probability measure v e being a QSD for P e , we have 

KS T d > f ) = e ~ Xet > v * > °- 

Thanks to (@J, \ € is decreasing to A = inf^ gC -°o(] 0j+oo [) (L o 0>0)^ when e goes to 0. As a 
consequence, 

P°(r a >t)= e- Xot , Vt > 0. 

In this case, the density of v with respect to dfi is an eigenfunction of L° with eigenvalue 
— A < 0, where (L )* is the adjoint operator of L° (this is a consequence of the spectral 
decomposition proved in 0, Theorem 3.2]). As defined, — A is at the bottom of the spectrum 
of (L )*. Thanks to @, Theorem 3.2], this eigenvalue is simple. Moreover, 0, Theorem 5.2] 
states that this QSD is equal to \im t -, +00 P°(A t G .\rg > t), what concludes the proof. □ 
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2.3 Diffusions with values in a bounded interval 

Theorem [1] is stated for 1-dimensional diffusions with values in ]0, + oof. However, most of 
the proofs can be easily adapted to diffusions with values in a bounded interval )a,b[, where 
— oo < a < b < +00, defined by the SDE 

dX t = dB t — q(X t )dt, X = x &)a,b[, 

and killed when it hits a or b. Here B is a standard Brownian motion and q G C 1 (]a,6[). 

More precisely, let us denote by P° the law of such a diffusion. For each e > 0, define P e 
as the law of a diffusion with values in }a + e,b — e[, driven by the SDE 

dX t = dB t — q(X t )dt, X = x e]a + e,b — e[, 



and killed when it hits a + e or b — e. As proved in [22j, there exists a unique QSD v e associated 
with P e . 

We define Q{x) = f^ a+b y 2 Q(y)dy and W(x) = q(x) 2 — q'(x). The counterpart of Theorem 
[1] under these settings is 

Theorem 9 Assume that the following hypotheses are fulfilled: 

Hypothesis 4 (HH1 ) W is uniformly bounded below by —C, where C is a positive con- 
stant. 

Hypothesis 5 (HH2) x t— > y^^ + c+i e ~ 2< ^^ or x ^ (x — a)e~®^ is integrable on a neigh- 
bourhood of a. 

Hypothesis 6 (HH3) x 1— > w ^ +c+1 e~ 2Q ^ or x 1— >■ (6 — x)e~ Q( - x ^ is integrable on a neigh- 
bourhood ofb. 

Then the family of QSD (i/ e ) is tight as family of measures on ]a,b[. Moreover, every limit 
point of the family (V e )o<e<i/2 ^ s a QSD for P°. 

Remark 5 Our aim isn't to develop this part, but we point out that to show that Proposition 
M remains valid, most of the arguments used to prove the key results 0, Theorem 3.2] and 
|3j, Theorem 5.2] can be adapted to these settings. 



3 Approximation of is e , QSD for P e 

We are interested in proving an approximation method for the QSD associated with F e . It 
will be sufficient to prove it for any diffusion (X t ) taking its values in ]0,1[ in place of ]e,l/e[, 
defined by the stochastic differential equation (SDE) 

dX t = dB t - q(X t )dt, X = x > 0, (16) 
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and killed when X t hits the boundary {0,1}. Here B is a real Brownian motion and q G 
C 1 ([0,1]). The law of X will be denoted by P. 

From |22j, the QSD of X is unique and equals the Yaglom limit. It will be denoted by v. 
For notational convenience, new notations have been defined for this section, which is totally 
independent of the previous one. 

Fix N > 2 and let us define formally the interacting particle process with N parti- 
cles described in the introduction. Let B 1 ,...,B N be N independent Brownian motions and 
(Xq,...,Xq) e]0,1[ n be the starting point of the process. 

• For each % G {l,...,iV}, the particle X % evolves in ]0,1[ and satisfies the SDE dX\ = dB\ — 
q{X\)dt (and then it is independent of the others) until r\ = inf{t > 0, X\_ = or 1}. 

• At time T\ = min{r 1 1 ,...,T^}, the path of a particle, denoted by i\ (it is unique), has a 
left limit equal to or 1. 

• A particle j\ is chosen in {1,...,N} \ {ii}. The particle i\ jumps on the position of the 
particle we set X^} := X^. 

• After time r 1; each particle X % evolves in ]0,1[ with respect to the SDE dX\ = dB\ — 
q(Xl)dt until t\ = mf{t > T\, X\_ = or 1}. At time Tl, all the particles are in ]0,1[, 
so that we have t\ > t\ for all i G {1,...,N} almost surely. 

• At time r 2 = min{r 2 1 ,...,r^ v } (which is then strictly bigger than ri), a unique particle %2 
has a path whose left limit is equal to or 1. 

• A particle j-i i s chosen in {1,...,A^} \ {i 2 }- The particle i% jumps on the position of the 
particle j 2 : we set X% := X-j?. 

• After time r 2 , the particles evolve independently from each other and so on. 

Following this way, we define the strictly increasing sequence of stopping times < T\ < 
t<i < r 3 < the time = linin^ooTn and the interacting particle system (Xl,...,Xf) for 
all t G [0,7-oof. The law of (X\...,X N ) will be denoted by F ipp . 

We can now state the main result of this section: 

Theorem 10 (X 1 ,...,X N ) is well defined, that means = +oo almost surely. It is geomet- 
rically ergodic, with unique stationary distribution M N . 

Let X N be the empirical stationary measure of the interacting particle process with N 
particles, that is the empirical measure of a random vector (a; 1 ,...,!^) G]0,1[ distributed with 
respect to the stationary measure M N of the process (X 1 ,...,X N ) . The sequence of random 
measures {X n )n>2 converges in law to the deterministic measure v, QSD of the process X. 

Subsection 13.11 is devoted to prove a coupling which ensures the non-degeneracy of the 
particles at the boundary. A consequence will be that r M = +oo almost surely. In Subsection 
I3.2[ the process is studied in finite time. We prove that the empirical measure of the process 
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(X 1 ,...,^^) at time t converges, when N goes to infinity, to the distribution of X t conditioned 
to not be killed until time t, which is P Mo (Xt G .\X t ^ d) (/io denotes the limit of the empirical 
measure of the process at time 0). We prove in Subsection 13.31 that the interacting particle 
system (X 1 ,...^^) is geometrically ergodic. We conclude by showing the convergence of the 
empirical stationary measure to the QSD v. 

3.1 Existence and non-degeneracy at the boundary 

One of the most important fact when studying the interacting particle system is that the 
particles don't degenerate at the boundary. This is an evidence if the particles are inde- 
pendent. In our case, we will prove a coupling between (X 1 ,...,X N ) and an other process 
(Y 1 ,...,Y N ) whose components are independent identically distributed and don't degenerate 
at the boundary and such that, for all % G {1,...,N}, 

< Y i < X 1 < 1 - Y i < 1 a.s. (17) 

With this construction, the process (X 1 ,...,^^) doesn't degenerate at the boundary, because 
each of its particles X % is contained in [Y l ,l — Y 1 }. This coupling will be useful in each step 
of the proof. 

3.1.1 Coupling's construction 

Define Q = swp xe pu\q(x)\ (we have Q < +oo by hypothesis) and fix i G {l,...,iV}. The 
process Y % is defined with values in [0,1/3] by the SDE 

dY'l = dWi - Qdt, F^O) = min {X%0),1 - JT(0),l/3}, 

with and 1/3 as reflecting boundary (see [5( for the definition of a reflected diffusion). The 
coupling inequality (JTTJ) is fulfilled at time t = 0. The Brownian motion W l will depend on 
B % and on the position of X % . 

If X\ belongs to [0,1/3], then X\ < 1/3 and the second part of the coupling inequality 
is satisfied, independently of the choice of W l . We only need to ensure that X\ stays bigger 
than F/. If X\ belongs to [2/3,1 [, we only need to ensure that it is smaller than 1 — Y t \ If it 
belongs to ]l/3,2/3[, the coupling inequality is obviously fulfilled, thanks to the reflection of 
Y* on 1/3. 

Assume that X\ is in ]0,l/3] and that the coupling inequality is fulfilled at time t. We 
have d(X\ - Y t l ) = dB\ - dW\ + (Q - q(X t ))dt, with Q - q(X t ) > 0. If we choose W i so 
that dB\ — dW\ = 0, then d(X l t — Y t l ) = (Q — q(X t ))dt is increasing with time almost surely 
and the coupling inequality remains fulfilled. In a similar way, if X\ is in [2/3,1 [, we have to 
choose W % so that dB\ + dW\ = 0. We will see in the proof of the coupling inequality that 
the jumps of X 1 and the reflexion of Y l on do not play any role in the coupling inequality 
(see Figure [2]). 

Let us define the Brownian motion W l in an explicit form. First, we build a sequence of 
strictly increasing times (^J n >o such that, for each n > 0, X\ e]0,2/3[ for all t G [^n^n+it 
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and X\ G] 1/3,1 [ for all t G [t\ n+1 ,t 

2n+2[- Define 
f = inf{te[0, + oo],X*e]0,l/3]}, 

t\ = M{te [t , + 00], XI e [2/3,1 [}, 

and, for n > 1, 

4 = inf{f G [4 n _ 1 , + oo],X*e]0,l/3]}, 
4 +1 = inf G [4, + 00], X* G]2/3,l]}. 

Depending on the position of X\ which is stated by the sequence (t n )n>2, we define W l by 

W; = -B\ for t G [0,f ], 

and, for all n > 0, 

w? = w t ; n + (s; - sjj for t g [w^i] 
= »Ci - - S L + J for * e ki>^- 

If linin^ootn < +00, then can be extended by continuity on [0, + 00 [, so that Y t l is well 
defined for all t G [0, + 00 [. 

The sequence (t l n ) ne ^ is a sequence of stopping times for the natural filtration of the 
process (X 1 ,...^ 1 *). Conditionally to the sequence (j n ) n >i (the particles which are chosen 
at each successive jump), (X 1 ,...^) only depends on (B 1 ,...,B N ). One can then apply the 
strong Markov property to (B 1 ,...^) at time t* n , Vi G {1,...,N} and Vn G N. 

As a direct consequence of the symmetry of a Brownian motion's law and of the strong 
Markov property applied to the iV-dimensional Brownian motion (B 1 ,...,B\...,B N ), W l is a 
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Brownian motion. Note that the sequence of stopping times depends clearly on the position 
of other particles and then on the B^s, for j G {1,...,N}. However, the Brownian motions 
{W*}j=i,...,N are independent processes. Indeed, the sequence (£„)„>o doesn't play any role 
in the law of W l and, to be convinced of that, one can compute the covariance matrix of the 
Gaussian variables W/, for any t G [0, + oo[, which is clearly a diagonal one. 

Proposition 11 The coupling inequality (ITT)) is fulfilled for allt G [0,roo[. 

Proof of Proposition [7JJ : Define the time ( = inf {0 < t < r^, Y t l > X\} and let us work 
conditionally to ( < r M , then we have, by right continuity of the two processes, Y£ > X]» a.s. 

We first show that ( is a jump time for the particle i. Assume the converse. If ( = 0, then 
X^ = Y£ and, if £ > 0, then X 1 and Y % are continuous in a neighbourhood [( — h£ + h] of (, 
where h > 0, and we have Y t l < X\, for all t G [( — h,([, so that, by left continuity Y^ < X^, 
and then Y£ = X^. Therefore, we have X£ g]0,1/3] and 3n > such that £ G [t2n^2n+i[> 
and, for h small enough, one can assume that ( + h < £2^+1- Then, for all t G [(,( + /i], 
— Y£) = (Q — q(Xl))dt + dL^ 3 , where Q — q{X l t ) > and L^ 3 is an increasing process 
due to the reflecting property of the boundary 1/3 for Y % (note that Y£ > 0, so that, for 
h > small enough, Y£ > for all t G [C,C + M)- Then X* — stays non- negative between 
times ( and ( + h, what contradicts the definition (. 

The time £ is then a jump time for the particle i. If X^_ = 0, then, by definition of £, 
= and, by left continuity of the process Y\ Y£ = 0, so that X^ = 0, what is impossible. 
Therefore X^_ = 1 > 1 — Y?, and, by existence of left limits for the two processes, 3t < ( 
such that X\ > 1 — Y t l . Define (' = inf {t > 0, 1 — Y t l < X l t }. We have then, conditionally to 

C < t-oo, c < c 

By symmetry, conditionally to the event (' < r M , we have ( < then ( < r M and 
therefore ( < Finally, ( = (' = r M almost surely. □ 

3.1.2 Existence of the interacting particle process 

Proposition 12 For all N > 2, the interacting particle system (X 1 ,...,X N ) is well defined, 
that is = +00 almost surely. 

Proof of Proposition [7H : Let N > 2 be the size of the interacting particle system and fix 
arbitrarily its starting point x e]Q,1[ n . We define the event C x = {tqo < +00}. 

Conditionally to C x , the total number of jumps is equal to +00. There is a finite number 
of particles, then at least one particle, say io, makes an infinite number of jump before t^. 
At each jump of io, a particle is uniformly chosen in {1,...,N}. By the law of large numbers, 
each particle is chosen infinitely often before r^. Assume that a particle, say jo, remains all 
the time in ]e,l — e[, with e > 0. io will jump on the position of jo infinitely often. Then it will 
come back from ]e,l — e[ to the boundary infinitely often in finite time, what is impossible. 
We deduce that, conditionally to C x , all particles of the interacting particle system are going 
near to the boundary, that is 

C x C 1 lim min (X l t ,l - X\) = ol , (18) 

t—tToo 
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for each particle i G {1,2,..., AT}. 

Using the coupling inequality of Proposition [TH we deduce from ffl8|) that 



C X C lim (Y t \...,Y t N ) = 




} 



Then, conditionally to C^, y 1 and F 2 are independent reflected diffusions with bounded 
drifts, which hit at the same time. This occurs for two independent reflected Brownian 
motions with probability 0, and then for Y l and Y 2 too, by the Girsanov's Theorem. That 
implies P{C X ) = 0. 

Finally, we have = +oo almost surely. 

Remark 6 One could hope to apply this method directly to the process with law P° studied 
in the first part of this paper. Unfortunately, it can be very difficult to show the existence of 
the process or the non-degeneracy at the boundary: the drift being unbounded, the law of 
the reflected diffusion used in this proof isn't absolutely continuous with respect to the law 
of the reflected Brownian motion and stay at all the time after hitting it. 

3.1.3 Non-degeneracy at the boundary 

For all r > 0, we define the open set D r =]r,l — r[. Let fi N (t,dx) (resp. fi' N (t,dx)) be the 
empirical measure of the system of particles (Xl)i = i t ___ t jsf (resp. (17)1=1,..., at), that is 



We will suppose that, at time 0, the sequence of empirical measures (fi N (0,dx))N> 2 satisfies 
the following non- degeneracy property, which ensures that the mass of fi N (0,dx) doesn't 
degenerate at the boundary, uniformly in N: 

Definition 1 The family of random probabilities {/j, N (dx)} is said to verify the non- degeneracy 
property if, for any e > 0, 



where D° =}0,r) U [1 — r,l[. 

From definition of (Vo)ie{i,...,iV}> the non-degeneracy of fi' N (0,dx) is the consequence of the 
non-degeneracy of [i N (0,dx). The end of the section is devoted to prove the following Propo- 
sition, which states that the non- degeneracy property is maintained over time: 

Proposition 13 Assume that (fi N (0,dx))N>2 satisfies the non-degeneracy property, then, for 
all T > and all e > 0, 



fi N (t,dx) = -zj^^fixiidx) and fi' N (t,dx) = — ^^^((ia;). 




lim lim sup P (^ N (D c r ) > e) = 



(19) 



lim lim sup P( sup n {t 1 D c r ) > e) = 0. 

JV^oo te[0,T] 
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Proof of Proposition [7^1 : Fix T > and e > 0. Because of the non- degeneracy of fjf N (0,dx), 
one can find a > such that 



(20) 



when N goes to oo. 

We want to apply a law of large numbers, but the Yq aren't independent and depend on 
the number of particles N. Let us define the diffusion Z % with values in [0,1/3], defined by 
the SDE 

dZ\ = dW? - Qdt, Z l (0) = a, 

with and 1/3 as reflecting boundary Here the W n are independent Brownian motions. The 
random processes l Z i &D c are independent, identically distributed with values in D([0,T],M). 
One can apply to them the following law of large numbers, proved in 24| : 



sup 

te[o,T] 



1 N 

i=l 



Prob 

N— »oo 



0. 



We have i£ (l^eD^) = -P {% % ^r), which tends uniformly to when r 
from it that 



lim lim P I sup — J2 1 z l eD ? > e/2 ) 



0. 



0. We deduce 



(21) 



For each number of particles N > 2, one can easily find a coupling between Z % and Y l 
where Z % < Yq implies Z\ < Y t l for all t G [0,T]. With such a coupling, we have 



1 A if 



i=l 



i=l 



By adding the contribution of the y* which starts in ]0,a[, we get 



^ n ^ TV 1 ^ 

i=l i=l i=l 



The limits ( 120]) and (ED) allow us to conclude the proof. 



3.2 Convergence in finite time 

Fix T > 0. This section is devoted to the proof of the following proposition, which states that 
the empirical measure process converges to the distribution of the process X conditioned to 
not be killed. 
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Proposition 14 Assume that (/jl (0,dx))Nm converges in law to the random probability 
measure n(0,dx) with respect to the weak topology and satisfies the non- degeneracy property. 

Then, VT > ; the measure processes ([i N (t,dx)) t <=[o,T] converge in law to (F(X t G dx\X t ^ 
9))te[o,T] in the Skorokhod space D([0,T],j\4i(]0,l[)) when N — > oo. Here «Mi(]0,l[) denotes 
the space of probability measures on ]0,1[ equipped with the weak topology. 

Proof of Proposition 14 ■' For all maps ip G Cn([0,11) vanishing on {0,1}, one can apply the 
Ito's formula to the semimartingale ip(X^) (see |2ll . Theorem 27.1]), whose number of jumps 
in [0,T] is finite almost surely : 

i>{x\) = 4>{x*) + jj\xl)dBi + £ U'(xl) q (xl) + -<(*:.)) ds 

+ E - ( 22 ) 

0<s<t 

Let us denote by (T^)i< n the increasing sequence of jump times of the particle i. We have 

E ^ X s) - = E MO - ^ X U (23) 

0<s<t 0<r^<t 

= E (2 4 ) 

0<ri<t 

because X\ G {0,1} and ip(0) = = 0. That implies 

J2^xi)-^xi)= e wy-^E^y) 

0<s<t 0<ri,<t \ j=l / 



(25) 



J = 1 0<T»<t 

By summing over i G {l,...,iV}, we obtain 



+ iV3I E (27) 



0<T„<i 

where Ai c (ip,t) is the continuous martingale 4 Xlili Jo ^ {X\_)dB\ and MPtyfi) is the pure 
jump martingale 



JV / TV 

m^*)=^e E ^y-F^iE^i- 

»=i o<rf<t V j=i 



(28) 
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Now, we interpret each jump as a killing. Then we introduce a loss of 1/N of the total 
mass at each jump: we look at the measure process fj, N decreased by a factor ^fei at each 
jump. More precisely, we set 

v N (t,dx) = [-^j » N (t,dx), 

where A N (t) = # Ui=i{ r nj — T n — t} denotes the total number of jumps before time t. 

Lemma 15 The sequence of measure processes (v N (.,dx)) N converges in law to P Mo (X. e dx) 
in the Skorokhod topology D({0,T],M(]0M)) ■ 

Proof of Lemma UM ■ Applying the Ito's formula to the semimartingale u N (t,i/j), we deduce 
from fl27|) that 

p N {t^) = v N {^)+ f\ N (s-^q + ^")ds + f S ~ 'dM c (^,s) 



o 



N 

+ (v N (Tn,i>)-V N (T n -,^)), 



0<T n <t 

Where we have 

s A N 

N — 1\ n 



. , ,'N-1\ K / A' - 1\ 

+ /i A, (r n -,^) 



with 



and 



N \ N 



A N , A N , s A N 



then 



N I V N N — 1 \ N 



, A N 



s A N 

(A^T n )-A^T n -)) 



iV-1 /JV- l x 



iV V iV 
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That implies 



+ 



2 

N — 1 



N 



E 

0<T n <t 



* fN -I 
N 



N - 1 

N 



dM c {ip,s) 

(M j (i/;,T n )-M j (<il>,T n -)) 



We deduce that for all smooth functions ty(t,x) vanishing at the boundary 



!/"(*,¥(*,.))- i/"(0,¥(0,.))- / 



9*(s,.) 9*(s,.) l<9 2 #(s,. 



9s 



+ 



<9x 



2 <9x 2 



-)ds 



N- 1 

N 



dM c (ty(s,.),s) 



+ 



N - 1 



A' 



E 

0<T„<t 



(^(*(r n ,.),r n ) - ^'($(r„-,.),r n -)) (29) 



Because (^r") As ~ < 1 &.s. and by the Doob's inequality , we have 

1 



E I sup 

te[o,T] 



iV 



djVt c (tf(t,.),s) 



< —Til — II 2 . 

" N ll dx llo ° 



(30) 



Note that the jumps of the martingale M? are smaller than j|||\l/||oo, then 



E 



E 

0<Tn<T 



N -I 

N 



2A T 



(A1 j (*(t„,.),t„) - X'(*(t„-,.),t„-)) 

E 







< 


-II* 




N 2 


< 


±11*11 




jyll II 



1^ 



.0<r„<T 



7\T _ i \ 2At «- 



iV 



2 



By the Doob's inequality, we have then 



TJ ( sup 

te[o,T] 



E 

0<T n <t 



A/" — 1 \ n 



A' 



(A^ J '(*(r n ,.),r n ) - A4 J '(*(T n -,.),T„-)) 



<^ll*llL (31) 



Define *&(s,x) = P t - S f(x), where / G C°°([0,1]) vanishes on {0,1}, and (P t ) is the semi- 
group associated with the diffusion X defined by (|T6|) . From Kolmogorov's equation (see 0, 
Proposition 1.5 p. 9]), 

d 1 

— W(s,x) + -AV(s,x) + q(x)W(s,x) = 0. 
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We deduce from (129]), (j30]) and (l3llh that 



E I sup 

te[o,T] 



v"(t,f)- I P t f(x)d^(0,x) 





< jjC(f), 



where C(f) is a positive constant, which only depends on /. For each map g G C°°([0,1]), 
one can set f r (x) = j r (x)g(x), with r > 0, where 7 r G C°°([0,1]) is equal to 1 on ]2r, 1 — 2r[ 
and vanishes on ]0,r[U]l — r,l[. Then 



^(t,*?) - I P t g(x)d^ N (0,x) 
Jo 



< 



v n (t,f r )-J P t f r (x)d^(0,x) 
+ \u N (t,(l- lr )g)\ 



P t {{l- lr )g) (x)d^{0,x) 



(32) 
(33) 



where ( 132]) (see Proposition [T3l) and (|33|) are going to when r tends to 0, uniformly in N, 
and f r G C°°([0,1]) vanishes on {0,1}. Then 



sup 

f€[0,T] 



u"(t,g)- I P t g(x)d^(0,x) 
o 



0. 



In particular, u (.,dx) converges in law to P^ (X. G dx). □ 
Let us conclude the proof of Proposition [TH From Lemma [151, 

(/(tjo.iD/^di))^^^ QP>,&oGx* e]o,i[),F^ dx) (x t g rfx)) ie[0 ,T] 

iv — >oo 

in the Skorokhod topology D([0,T),R x Mi()0,l[)). That means 
'iV- 1 



N 



v N (t,dx) 



*e[o,T] 



law 



P M (0,^)(^ e]0,lD: P M(0,«&:)(^ e ^)) te[ 0,T]- 



The process P^o.dx) (^"- s]0,l[) never vanishing almost surely and the limit process being 
continuous almost surely, we have 



N-l 

N 



u N (t,dx) 



t€[0,T] 



law 



te[o,T] 



in the Skorokhod topology £)([0,T],A^i(]0,l[)). That means 



(/i W (t,rfx)) ie[0 ,Ti (P M (o,<fa)(-X't ^ dx|X t ^ d))te[ 0t T\ 



law 



N^oo 



in the Skorokhod topology D([0,T],A / li(]0,l[)). The proof of Proposition [LH is then complete. 
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3.3 Existence and convergence of the empirical stationary mea- 
sures X 

For each N > 2, we say that the interacting particle process (X 1 ,...,^) is exponentially 
ergodic, if there exists a probability measure M N on ]0,1[ such that, 

\\P x ((X},...,X t N ) G .) - M N \\ TV < C{x)p\ Vx Vt G R+, (34) 

where C(x) is finite, p < 1 and ||.||rv is the total variation norm. In particular, M N is a 
stationary measure for the process (X 1 ,...^). When M N exists, we denote by X N the 
empirical stationary measure associated with M N , that is a random probability which is 
distributed as jr^2 i=1 S Xi , where (x 1 , ...,x N ) is a random vector in jO,!^ distributed with 
respect to M N . 

In a first time, we prove that for all N > 2, the interacting particle process with N 
particles (X 1 ,...,^^) associated with X is exponentially ergodic. We conclude by proving 
that (X n )n converges in law to the unique QSD of X. 

3.3.1 Exponential ergodicity 

Here N > 2 is fixed. We are interested in proving the following result, which is the first part 
of Theorem [10] 

Proposition 16 The interacting particle process (X 1 ,...,^) with law P w is exponentially 
ergodic. 

Proof of Proposition [Tb] : We focus on the 1-skeleton of the interacting particle process with 
N particles, which is the Markov chain (X^,...,X^) n£ ^. Thanks to [a, Theorem 5.3 p. 1681], 
exponential ergodicity of (X 1 ,...,X N ) will be obtained as soon as the associated 1-skeleton 
is geometrically ergodic, which means that it exists a probability measure n N on jOT^ such 
that 

\\P a {(Xl...,X?) G .) — n N \\ TV < C (x)p n Ql Vx G]0,ir, Vn G N, 

where Cq(x) is finite and po < 1- 

To prove the geometrical ergodicity of the 1-skeleton, let us introduce the following defi- 
nition: 

Definition 2 C cjOT^ is said to be a small set for the Markov chain (X^,...,X^) neN 
if, for some nontrivial probability measure $ and some n > 1, e > 0, the n-step transition 
probability kernel P n (x,A) := P*? P ((X^,...,X^) G A) satisfies, for all xeC, 

P n (x,A) >e&(A), AeB(]0,l[ N ). 

Lemma 17 All compact set C = [r,l — r] N (with r > 0) is a small set for the 1-skeleton. 
Moreover, Ek > so that 

supE x (ht T c) < oo, (35) 

zee 

where r' n is the return time to C . 
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Proof of Lemma 17: Fix r > and let T be the event "the process (X 1 ,...,X ) has no 



jumps between times and 1". Define p = inf, rg [ r . il _ r ,]iv P*f p (jF). Thanks to the coupling with 
(Y 1 ,...,Y N ), we have p > 0. Conditionally to the event J 7 , the particles of the interacting 
particle process are independent from each other. 

Let us study d l (dx) = P^ x ^(Xl G dx and J 7 ). The law of X 1 conditionally to T 
is the same as the law of X conditioned to not jump, because, given this last event, the 
process X 1 doesn't depend on the other particles. Thus the probability of l X\ G dx and J 7 " 
is P-j^Xi G dx). The law of X\ has a density p\(xi,y) with respect to the Lebesgue's measure 
and pi(xi,y) depends continuously on x\ and y. It only vanishes when y = or 1. Then 

inf Pi(xi,y) > 

(x 1 ,y)e[r,l-r]x[r,l-r] 

Denoting this minimum by e', we have, for all X\ G [r,l — r], d l (dx) > e'l[ r ^ r ](x)dx. 
Conditionally to the particles are independent from each other, so that 

N 

where P!^ x JV )(^i e ^l^ 7 ) is greater than P^ Xn ^(X{ G dyi and JF) and then greater 
than d l (dyi). Finally, we have 

P (xL..,xjv)^ X i 1 '"-' X i V ) G d yi- d VN\^) > e ,N l [r!l - r] N(y 1 ,...,y N )dy l ...dy N , 

so that [r,l — r] N is a small set. 

For all x &]0,1[ N , and all n > 1, the probability of being in C at time n + 1 starting from 
a; at time n is bounded below by the probability pc > for (Y 1 ,...,Y N ) to enter C at time 
n + 1, starting from at time n. Hence, at each time n > 1, (X 1 ,...^^) returns to C at time 
n + l with a probability greater than pc > 0. That implies that the return time to C for the 
1-skeleton of the interacting particle process with X particles is bounded above by a time 
of geometrical law, independent of the starting point x G]0,1[ , and then satisfies condition 
(USD. □ 



The chain (X^,...,X^) neN is aperiodic. Moreover, if the Lebesgue measure of a subset 
A cjOjlf"^ is strictly positive, then F 1 p p {ta < oo) > for all x e]0,1[ n , where ta is the first 
hitting time on A for the chain (X^,...,X^) ng N- Thanks to 0, Theorem 2.1 p. 1673], if such 
a Markov chain has a small set which satisfies (|35j) . then it is geometrically ergodic. As a 
consequence, Lemma [17] allows us to conclude the proof of Proposition [T6j 



3.3.2 Convergence to the QSD 

We are interested in proving the following result, which is the second part of Theorem [TU1 

Proposition 18 The sequence of random measures {X N ) N>2 converges in law to the deter- 
ministic measure v, QSD of the process X. 
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Proof of Proposition HE ■' For each r g]0,1/4[, we define 7 r as a non- negative continuous 
bounded function from ]0,1[ to R, equal to 1 on D 2r and equal to on D c r . We have 

X N (D c r )<X N (l- lr ),Wre}0,l/4[. 

Thanks to Proposition [T6], the sequence of random measures (fi N (t,dx))N>2 converges in law 
to X when t tends to +00. That implies 

E (n N (t,l - lr )) — E (X N (1 - Tr )) , Vr G]0,l/4[. 

We denote by /j,' N (t,dx) the empirical measure of (Y t l ,...,Y t N ). Let us choose 7 r monotone on 
]r,2r[ and ]1 — 2r,l — r[. From the coupling inequality, 

» N (t,l- lr )<»' N (t,l- lr ) 

for all t G M + and r > 0. Then 

E(» N (t,l-7r))<E(n' N (t,l- lr )), 

which tends to when t tends to +00 and r to 0, uniformly in N. As a consequence, 

E{X N ){l- lr )— >0, 

r— >0 

where E(X N ) is the deterministic measure defined by £ , (Af JV )(^4) = £'(A' Ar (y4)), for all mea- 
surable set A. That yields 

E(X N )(D*)— >0, 

r— >0 

uniformly in AT. The family of intensity measures (E(X n ))n is then tight. This is a sufficient 
condition for the family of random variables (X ) to be tight, as shown in 16|, Corollary 
2.2]. We conclude that it exists a sub-sequence (A^^) which converges in law to a random 
probability measure X . 

Choose /j, N (0,dx) = X (dx). The non-degeneracy property is fulfilled. Thanks to Propo- 
sition dn 

(^ N \t,dx)) te[0 ,T] ^ (fx(x t e dx\x t ^ d)) te[0 , T] , vr > 

in the Skorokhod topology D([0,T],A4i(]0,l[)). The limiting process (Fx(X t G dx\X t ^ 
d)) t £[o,T] being almost surely continuous, 

fi m \t,dx) ^ F x {X t G dx\X t ^ d), Vt > 0, 

N~^oo 

with respect to the weak topology of A4i(]0,l[). By stationarity, the random probability 
measures ' (t,dx) and X*W have the same law. Making N tend to 00, we deduce that 
^x{Xt G dx\X t 7^ d) and X have the same law too. This looks like a QSD, but X is a priori 
a random measure and we need the following result to conclude. 

Lemma 19 For all m G A4i(]0,l[) ; 

lim P m {X t G A\r >t) = v{A), 

t^oo 

where r is the killing time of the process X and v its unique QSD. 
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Proof of LemmaUM ■' Let m be a probability measure on ]0,1[. 3Ao > 0, </>o and </> two con- 
tinuous maps vanishing on and 1 such that, for all x G]0,1[ (see R.G. Pinsky's explanations 



22l . Hypotheses 2 and 3]): 

lim e Xot P x (r > t) = Ci<f) (x), (36) 

lim e Xot P x (X t eA,r>t) = C 2 o (x) / My)dy. (37) 

t-^oo J A 



Here e Xot P x (r > t) is uniformly bounded above in the variables t and x (see [3|, Proof of the 
equality 7.2, p27]), then, by dominated convergence, one can integrate with respect to m 
under the limit in ([36]), 



lim e Xot P m (r > t) = d I 

t^oo J D 



4>o(x)m(dx) 



The same holds for fl37l): 



lim e Xot P m (X t eA,r>t) = C 2 I (j) ( x ) I Mv)dy m(dx 



t^oo 



D J A 



Then, by Fubini's Theorem, 



Pm{X t EA,T>t) C 2 

^ P rn {r >t) = C[ I My)dy > 



A 



that is, from flUD and (07]) with A =]0,1[, 

P m (X t € A, r>t) = f A My)dy 
*- +D0 Pm(r>t) J D My)dy 

which is nothing else but u(A) (see [221, Proposition 1.10]). □ 

Thanks to Lemma [T9| Fx{X t G dx|X t ^ d) converges almost surely to the Yaglom limit 
when t — ► +oo, and so do Af. Finally, A? is the unique QSD of the process and the proof is 
complete. □ 



4 Numerical applications 

4.1 The logistic case 

We apply our result to the logistic Feller diffusion with values in ]0, + oo[, driven by the 
stochastic differential equation 

dZ t = \[Z t dB t + (rZ t - cZ?)dt, Z = z > 0, (38) 
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and killed when it hits 0. Here B is a 1- dimensional Brownian motion and r,c are two positive 
constants. 

We define P° as the law of 1\fZ., which is killed at and satisfies the SDE 



dX t = dB t 



2X t 



+ — ± ] dt, X Q = x e]0, + oof. 



For each e > 0, we define the law P e and denote its QSD by v e . 

As proved in (HI) and (H3) are fulfilled in this case. Thanks to Theorem [T] and 
denoting by v the Yaglom limit associated with P°, we have 



law 
U e >V. 



In the numerical simulations below, we set e equal to 0.001. 
By Theorem [TOj we have 

X N — > v f , 



where X N is the empirical measure of the system studied in Section [3j In the numerical 
simulations, we set N = 1000 and, because of the randomness of X N , we approximate 
E{X N ) using the Ergodic theorem: we compute X^=i°° ^N(t,dx). The graphic below 



(see Figure [3]) shows this approximation for different values of r and c. 

As it could be wanted for, greater is c, closer is the support of the QSD to 0. We thus 
numerically describe the impact of the linear and quadratic terms on the QSD. 



"i 1 — 

r=l,c=10 

r=l,c=l 

r=10, c=l 




012345678 
Figure 3: E(X N ) for the diffusion fl38|) . with different values of r and c 
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4.2 The Wright-Fisher case 

We illustrate the result of Section 12.31 by an application to the Wright-Fisher diffusion with 
values in ]0,1[ conditioned to be killed at 0. This diffusion is driven by the SDE 



dZ t 



%(1 - Z t )dB t - Z t dt, Z = z e]0,l[, 



and killed when it hits (f is never reached). In 13], the author proves that the QSD of this 
process exists and has the density 2 — 2x with respect to the Lebesgue measure. 

Define P° as the law of X = arccos(l — 2Z ), where Z is defined as above. P° is the law 
of the diffusion with values in ]0,ir[, driven by the SDE 

dX t = dB t - 1 ~ 2C °l f dt, X = x g]0,tt[, 
2 sin X t 

killed when it hits (it is never reached). For all e e]0,7r/2[, define P e and v e as in Section 

m 

The drift of the diffusion is q(x) = 1 ^si^xf* , Vx G]0,7r[. Let us show that it satisfies the 
hypotheses of Theorem [91 
We have, Va; G]0,7r[, 



q{x) 



q'(x) 



8 cos 2 x + 2 cos x + 2 cos x + 1 
4 sin 2 x 



+ 1 



which is positive and tends to +oo both in 0+ and 7r— . It implies that hypothesis (HH1) is 
fulfilled. 

For all x e]0,7r[, 



Q(x) 



q{y)dy 



tt/2 

it 

2 



X 

tan- 



+ In Isinxl 



In 2 sin — 

2 



„ 3 
X\ 2 



1 

X\ 2 



cos- 



That implies 



then 



and 



-Q(x) 



2 (sin a)' (cos f) 



xe 



-Q{x) 



7T — x e 



-Q(x) 



7T — X 



Finally, hypotheses (HH2) and (HH3) are satisfied and Theorem [9] can be applied. 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
Figure 4: Wright Fisher case, e = 0.001 and iV = 1000 



In the following numerical simulation (see Figure HJ), we set e = 0.001 and iV = 1000. We 
compute E(X ), which is an approximation of v e , and then of z/, with the method used in 
the logistic case (see Part 14.11) . 

The simulation is very close to the QSD (2 — 2x)dx, which shows the efficiency of the 
method. 

Remark 7 In the simulations, we have chosen to simulate a system with N = 1000 particles, 
because of empirical consideration. The question of convergence speed will be studied in 
further works. 
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